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5 ■ Abstract 

-)— > 

I . Variational Monte Carlo method is used to calculate ground state properties 

^ . of fe droplets, containing 70, 112, 168, 240, 330, and 728 particles. The 

■ 

resulting particle and kinetic energy densities are used as an input in the 

^ . Feynman-Lekner theory for ^e impurities. The kinetic energy density of '^He 

atoms and the energy of the ^e surface states are compared with the results 

of previous phenomenological calculations. 
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I. INTRODUCTION 

Helium droplets have attracted some interest in recent years. A major motivation is the 
fact that they represent a prototype of finite size quantum systems, behaving as fluid of 
strongly interacting structureless particles. Several theoretical schemes have been developed 
for pure helium droplets [0-|^. The comparison between theory and experiments is still 
elusive, mainly because helium droplets are so weakly bound objects that their experimental 
characterization is very difficult. An interesting approach, from this point of view, is the use 
of atomic and molecular impurities as probes 0J§]- Theoretical calculations for droplets with 
impurities have been recently done with both Monte Carlo P| and density functional methods 
|TO| , [TT|| . Impurities heavier than helium atoms are expected to have small zero point motion. 



so that they can be treated as classical objects in a quantum fluid [|TT|. This is certainly 
not true for ^e and hydrogen impurities, for which a full quantum mechanical treatment is 
needed. The case of "^e impurities is particularly interesting from the theoretical viewpoint. 
Since ^e and ^e atoms interact through the same potential, the properties of their mixed 
systems are determined only by quantum effects, i.e., the different statistics and the different 
zero point motion. One '^e atom, being lighter than ^e , tends to move in regions of low 
^e density. This is the origin of the so called Andreev state of ^e on a ^e liquid-vapour 
interface at low temperature |]T^|-|I^. The same behavior is observed in ^e films on solid 



substrates |T8|, where the layer structure of the "fte density produces a rich variety of ^e 



states ||l9|- p^] . Predictions on the Andreev states of ^e atoms on ^e droplets have been 



already given |]10| using a phenomenological density functional, as well as the variational 
Feynman-Lekner theory. A key ingredient in the latter approach is the ^e kinetic energy 
density, for which a simple approximate expression was proposed in Ref. [|T5| . 

In the present work we calculate ground state properties of pure ^e droplets by means 
of a variational Monte Carlo ( VMC) method. We follow the same procedure as in Ref. p| , 
but with a different parametrization for the variational wave function. A detailed discussion 
about the results for the ground state of the droplets is given elsewhere p3| . Here we present 



the first microscopic results for tlie kinetic energy density, wliicli are sliown to be in good 
agreement witli tfie predictions of Ref. |T^. Finally, we use tlie ground state properties of 
^e droplets as an input in the Feynman-Lekner theory for ^e impurities. We calculate the 
binding energy of the ^e surface states for several ^e droplets. The results are extrapolated 
to estimate the binding energy of ^e atoms on a planar ^e surface, in good agreement with 



the experimental value [|14 



The work is organized as follows. Section II contains the short description of the VMC 
method and the results for the kinetic energy density of ^e droplets. The Feynman-Lekner 
approach is briefly introduced in Section III, where the results for the Andreev states are 
also discussed. Section IV is devoted to conclusions. 

II. PARTICLE AND KINETIC ENERGY DENSITY OF PURE ^He DROPLETS 

First, we want to calculate the ground state energy, density profile, and kinetic energy 
density for droplets of given number of particle A^. We use the VMC method. In the 
variational approach a suitable form for the many body wave function \l/t,(ri, . . . , tat) is 
choosen, containing a set of parameters that are varied in order to minimize the energy of 
the system. The Hamiltonian has the usual form 



where m^ is the mass of ^e atoms. Aziz HFDIIE2 potential p^ is used as the pair inter- 
atomic potential vijij). 

Several forms of the variational wave function of the pure ^e liquid have been proposed 
in the literature [ll|,^,[5|J^j23|l . They are expected to give very similar results in the context 



of the present work. We use the formalism of Ref. p3| , where the wave function is taken of 
the same form as in Ref. 0, 

^v=X{h{ri) n h{rij) n /3(ni,rjfc,rfci) , (2) 

i<N i<j<N i<j<k<N 



but with an improved version of the correlation functions /i, /2, and /s. This choice ensures 
a correct asymptotic behavior of the wave function far outside the droplet, which has an 
important role for the ^e and hydrogen impurity states. 

Once the optimal wave function is found, it is used to calculate the one-body density 
and the kinetic energy density, given by 



p(ri) = N J dr2...drN\'i^, 



(3) 



and 



A^ 



nri) 



dv2... TAr^* 



-h' 



-V\ ^, 



(4) 



p(ri) J "'^ " " \2m4, 
respectively. 

We consider droplets with 70, 112, 168, 240, 330, and 728 particles. The resulting density 
profiles and the kinetic energy density are shown in Fig. |l|. It is worth noticing that both 
Monte Carlo ||l|J^J^-|6|p3| and density functional P,pT| calculations predict an extrapolated 
free surface thickness (distance between the points where the density is 90% and 10% of the 
bulk value) of the order of 6 ^ 8 A, which is consistent with our results. 

The kinetic energy density is shown in the lower part of Fig. |1[ It decreases smoothly 
from the inner value, close to the value of the kinetic energy per particle in bulk liquid, to 
the asymptotic limit ^^{N), i.e., the chemical potential of the drop containing A^ ^e atoms. 
This limit follows from the behavior of the one-body factor which in our parametrization 
121 is 



/i(r -^ oo) oc y p{r -^ oo) oc -exp 



2m4|y[X4| 



(5) 



which dominates the r — >■ cxd limit of the wave function, yielding 



T[r — * ooj = /i4 



(6) 



In Ref. |]T5[ the approximate expression 



T = To 



.Po. 



2m4 ^ 



(7) 



was proposed for r, where tq and po are the ground state kinetic energy per particle and 
the particle density in bulk ^He at zero pressure, while n is a phenomenological parameter. 
Equation (|^) is an interpolation between the expected behavior of r in the two opposite 
limits p —^ po and p —>■ 0. It has never been checked so far with microscopic calculation. In 
Ref. []T3[ the values of tq and n were fixed in a phenomenological way to reproduce known 



properties of '^e impurities in bulk liquid "fte , while the saturation density was taken 
directly from experiments (tq = 13.34 K, n = 1.76, po = 0.365cr~^, where a = 2.556 A). 
The same quantities can be calculated microscopically. For instance, the kinetic energy per 



particle in bulk liquid at several densities was calculated in Ref. [^ using the variational 
wave function (0). From those data one extracts po = 0.362o"~^, tq = 14.52 K and n = 1.77. 
In Fig. 1^ we compare the VMC results for r, for droplets with 112 and 728 atoms, with the 
ones of the approximated formula (^ . The solid line corresponds to the variational values of 
To, n and po, while the dashed line corresponds to the phenomenological parameters of Ref. 
p!5| . In both cases we have used Eq. (|^ with the VMC density of the corresponding droplets, 
and the high frequency statistical fluctuations have been filtered out in the calculation of 
the second derivative. The figure reveals that the approximate formula (J^) works very well, 
especially for large droplets. 

III. LEKNER-FEYNMAN THEORY FOR ^He IMPURITIES 

Once the ground state properties of pure ^e droplets are obtained, one can calculate 
the energy and wave function of one ^e impurity. Consider a droplet consisting of (A^ — 1) 
^e atoms and one ^e impurity atom. Since the interatomic potential is the same for ^e 
and ^e , one can write the Hamiltonian in the form 

H = Ho + Hj , (8) 

where Hq is given in Eq. (|1|) and 
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A realistic trial wave function would be 

N 

v|/ = /(ri)n^(ri,r,,)vl/o , (10) 

1=2 

where \&o is the ground state of the Hamiltonian Hq, and / and F are variational functions 
to be determined. In principle, the optimal / and F can be found by minimizing the total 
energy of the system (\l'|if |\1'), which can be done either by solving resulting Euler equations 
or by direct Monte Carlo simulations. 

A significant simplification is obtained by taking F{ri,ri) = 1 [^. This corresponds to 
the assumption that the correlations between '^e impurity and ^e atoms are the same as 
those among ^e atoms. The Euler equation for / then becomes 

- T^Vi ■ (p(ri)Vi/(rO) + f ^ - l) /(rOr(rOMri) - A/(ri)p(ri) = , (11) 

where p and r are the particle and kinetic energy densities of the pure '^e droplet (see Eqs. 
(0) and (^)). The same equation can be rewritten in the form of a Schrodinger equation 

-^V\ir) + Vsiv)x{r) = exir) , (12) 

21713 

where 



X(r) = /(r)^p(r) (13) 

is the impurity wave function, and 



^3W^(!=l-i).W + ^^ + ,, (14) 

Vmg J 2m3 /p(r) 



is an effective potential seen by the ^e atom. The chemical potential /i4 of the droplet of 
the considered size is added for convenience, so that the eigenvalue e in Eq. (|T2|) is referred 
to the state in vacuum. In principle, an additional term should appear in the effective 
potential V3(r) because the exact ground state of pure ^e droplet \I'o is approximated by 
the variational wave function \l/„. This term is expected to be small due to the closeness of 
these two wave functions, and is neglected in the present treatment. 



The assumption F(ri,rj) = 1 in Eq. ([I0[) corresponds to the Feynman-Lekner theory 
, aheady used in the past to predict the properties of the Andreev state on the ^e free 
surface p!3|Jl5|] , on films ||22|, as well as on droplets [|T0|. Unlike in the previous calculations, 



where phenomenological p and r were used, here we take them from ah initio calculations. 
In order to find the energy and the wave function of the impurity on a droplet Eq. (|T2[) 
is rewritten in spherical coordinates: 

-yXnl +\Vz + 7, o Xnl = ^nlXnl ■ (15) 



2777-3 dl^^ V ^"^3 '"^ 

The crucial point is that the effective potential V3 has a well on the surface of the droplet, 
so that the lowest eigenstates of the Schrodinger equation are localized on the surface. The 
potential well originates from a balance between the excess kinetic energy of an ^e atom in 
the bulk with respect to the one of ^e atoms, which tends to push the ^e atom out, and 
the He-He interaction, which binds the ^e atom to the liquid. A typical situation is shown 
in Fig. ^ for a droplet of 112 atoms. Results for the energy of the lowest {n = 0) impurity 
states on six droplets are given in Fig. ^ as a function of A^~^/^. We note that a suitable 
smoothing procedure has been apphed to the VMC data in order to avoid spurious effects of 
statistical fluctuations of p in the calculation of the second derivative. This affects mainly the 
external tail of V3 where the density is very small. The error bars in Fig. ^ correspond to the 
consequent inaccuracy in the results, estimated by choosing different smoothing methods. 
The extrapolation to the case of the planar surface (A^ — > cxd) can be done by a linear fit, 
even if the accuracy of the fit is relatively poor. We obtain e(oo) ~ —4.9 K, rather close to 
the experimental estimate e(oo) = (—5.02 ± 0.03) K quoted in Ref. fl^. 



In Fig. ^ the results of the VMC calculation are compared with predictions of Ref. 
||T0| (full line) obtained using the Feynman-Lekner theory with approximation (0) for r and 
density profiles of Ref. 0. The significant difference between the two predictions is mainly 
due to the different density profiles. The curvature of the surface profiles in the outer region 
is underestimated in the density functional calculations of Ref. [Q with respect to the VMC 
results (see Fig. H). This makes the potential well for the impurity wider, and the binding 



energy lower. The density profiles calculated with a more recent density functional p7|JTl 
are sharper than the ones of Ref. |0 and the corresponding predictions for the ^e binding 
energy are closer to the VMC results (dashed lines in Figs. ^ and |^). Part of the remaining 
discrepancy is due to the use of the approximated expression (|^ for r; we have checked that 
the corresponding effect on the ^e binding energy is small (about 0.1 K). 

By solving the Schrodinger equation (|l^) one finds, above the lowest eigenstate, a spec- 
trum of states with different principal quantum number n and angular momentum /. The 
general features of the spectrum are the same as in Ref. |]lOl- In the limit of an infinite 



droplet the states with different / coincide with those of a 2-dimensional Fermi gas. 

To conclude, we stress again the idea of the Feynman-Lekner approach and its limits. 
The idea is that the two-body correlations between the '^e impurity and ^e atoms are 
taken to be the same as the ones between the ^e atoms. The ^e wave function is then 
expressed by means of the factor /(r) in the many body wave function (plQl). The form of 
/ is derived variationally, by solving a Schrodinger-like equation. The lowest eigenvalues 
are localized on the surface of the droplets, as an effect of the different mass of ^e and 
^e . To improve the Feynman-Lekner theory one should account for the fact that ^e -^e 
correlations differ from the ^e -'fte ones, i.e., one should take -F(ri, Fj) 7^ 1. One possibility 
is to choose F = (1 + 0), where (p is small. Keeping the terms in the expansion of the 
total energy up to the second order in is equivalent to taking into account the energy 
of an elastic deformation of the residual droplet due to the presence of the impurity. Such 
corrections have already been studied in Refs. [Q and |2^ in the case of liquid "fte films. 
The problem of ^e -^e correlations, in the context of variational calculations, has been also 
discussed in Ref. P^JSO] for bulk ^e -^e mixtures. Work in this direction is in progress. 



IV. CONCLUSIONS 



We present the variational Monte Carlo (VMC) calculation for ^e droplets with and 
without ^e impurities. We have used the variational wave function of Ref. M modified as 
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in Ref . p3[ . The results for the density and energy of six droplets, from A^ = 70 to A^ = 728, 
are close to the predictions of previous calculations. We have calculated the kinetic energy 
density, for which only approximate estimates have been given so far. Our results show that 
the analytic formula given in Ref. [^ for the kinetic energy density (Eq. (7)) works very 
well, and can be safely used in calculations involving the "^e surface. 

In the second part of the work we apply the Feynman-Lekner approach to study ^e 
impurity states on ^e droplets. The impurity wave function turns out to be localized on the 
surface of the droplets, as expected. The lowest value of the binding energy is almost linearly 
dependent on A^~^/^ and the value extrapolated to A^ — > cxd is close to the experimental 
binding energy of the Andreev state on a planar surface. We have also compared our 
results with the predictions for the '^e binding energies obtained using phenomenological 
calculations of the density profiles and the kinetic energy density. The general trend of these 
earlier results is similar to the one of VMC calculations. 
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FIGURES 
FIG. 1. Particle density and kinetic energy density for droplets of 70 (squares), 112 (plus 

signs), 168 (circles), 240 (crosses), 330 (triangles), and 728 (diamonds) atoms. 

FIG. 2. Kinetic energy density for droplets of 728 and 112 atoms. Circles: VMC; solid lines: 
Eq.(7) with po = 0.362cj^^, tq = 14.52 K and n = 1.77; dashed lines: the same equation with 
Pq = 0.365CJ-3, To = 13.34 K and n = 1.76. 

FIG. 3. Lowest eigenstate of one ^e impurity on a droplet of 112 atoms. Solid line: ^e wave 
function; squares: '^e density; circles: effective potential V3. 

FIG. 4. Binding energy of the lowest impurity state as a function of N"^'"^. Points with error 
bars: Feynman-Lekner theory with r and p from VMC calculations; solid line: Feynman-Lekner 
results of Ref. [10]; dashed line: same theory but with more recent density functional calculations 
for the density profiles, as in Ref. [11]. 

FIG. 5. Surface profile for a droplet of 112 atoms. Points: VMC results; solid line: density 
functional results of Ref. [3]; dashed line: density functional results of Ref. [11]. 
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